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ABSTRACT 

We report the discovery by the CoRoT space mission of a new giant planet, CoRoT-20b. The planet has a mass of 4.24 ± 0.23 M Sup and a radius 
of 0.84 ± 0.04 Rj up . With a mean density of 8.87 ± 1.10 gem -3 , it is among the most compact planets known so far. Evolution models for the 
planet suggest a mass of heavy elements of the order of 800 M s if embedded in a central core, requiring a revision either of the planet formation 
models or of planet evolution and structure models. We note however that smaller amounts of heavy elements are expected from more realistic 
models in which they are mixed throughout the envelope. The planet orbits a G-type star with an orbital period of 9.24 days and an eccentricity of 
0.56. The star's projected rotational velocity is v sin i = 4.5 ± 1 .0 km s , corresponding to a spin period of 1 1 .5 ± 3. 1 days if its axis of rotation is 
perpendicular to the orbital plane. In the framework of Darwinian theories and neglecting stellar magnetic breaking, we calculate the tidal evolution 
of the system and show that CoRoT-20b is presently one of the very few Darwin-stable planets that is evolving towards a triple synchronous state 
with equality of the orbital, planetary and stellar spin periods. 



Key words, stars: planetary systems - stars: fundamental parameters - techniques: photometry • 
troscopy 
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1. Introduction 

The existence of a planet population at very short orbital dis- 
tance, a < 0.1AU typically, with its wide range of orbital and 
physical properties is an intriguing phenomenon. In-situ forma- 
tion of such massive bodies so close to their host-star at a lo- 
cation where the circumstellar material is depleted and warm, 
appears indeed highly unlikely. Planet migration from further 
away distances where solid material is abundant, triggered by 
gravitational interactions, is invoked to account for this pop- 
ulation. The exact process is still under debate but two ma- 
jor mechanisms are put forward: gradual plane t migration due 
to interaction wi t h the circumstellar gas disk dLin et al.L Il996t 
Papaloizou et al., 2007) or planet-planet or pla net-companion 
star interactions combined with tidal dissipation (iRasio & Fordl 
1996; iFabrvckv & Tremaind l2007t iNagasawa et aUl2008l e.g"7). 
Whatever the exact nature of the formation path, these planets 
have undergone a significant orbital evolution since the time of 
their formation. Their current properties provide valuable hints 
helping to better understand their orbital evolution, especially 
the planet's orbit eccentricity and spin-orbit alignment of the 
system. 



* The CoRoT space mission, launched on December 27th 2006, 
has been developed and is operated by CNES, with the contribution 
of Austria, Belgium, Brazil , ESA (RSSD and Science Programme), 
Germany and Spain. 



Up to now, transit surveys have been more sensitive to 
planets at very short orbital period. The situation has recently 
started to change thanks to the extended temporal coverage of 
ground-b ased transit surveys and the adven t of s pace missions, 
CoRoT dBaglin et al l 120091: iDeleuil et all 120111) and Kepler 
(Boru ckietalll2010h . As a consequence, the number of plan- 
ets with orbital periods greater than a few days has significantly 
increased over the past two years. While the mean eccentric- 
ity for the close-in planets is close to zero, the transiting giant 
planets at larger orbital distance display a much wider range in 
eccentricity, a picture more consistent with the sample of plan- 
ets found by radial velocity surveys. These trends favor a third 
body induced migration with tidal circulari zation of an initial 
eccentric and po s sibly high-obliquity orbit (W inn et all [2010; 
Matsumu ra et all 1201 PL IPont et all 1201 ll e.g.). A consequence 
of this orbital evolution is the tidal destruction of the planet 
which spirals down onto the star in th e life time of the system 
dGonzalezl Il997t iJackson et all [2009), a dramatic destruction 
that app ears being the fate of th e vast majority of the transiting 
planets dMatsumura et alll2010l) . 

In this paper, we report on the discovery of CoRoT-20b, a 
new member of the hot-Jupiter class population. The planet tran- 
sits its G-type parent star every 9.24 days, along an orbit with 
a high eccentricity. The CoRoT observations are presented in 
Section [2] The accompanying follow-up observations and then- 
results are described in Section|3]and Section[5]for the host-star 
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Fig. 1. The 24.278-days long CoRoT-20 reduced light curve at a constant 51 2-sec time sampling. 



analysis. The final system parameters are derived in Section [4] 
We then discuss the properties of this new planet in Section [6] 
We investigate its orbital evolution and fate but also its internal 
structure that rises new questions on the nature of such compact 
object. 



2. CoRoT-20b Light curve 

The planet has been discovered in one of the fields observed by 
the C0R0T satellite (iBaglin et all 12009b iDeleuil et all 1201 ll) in 
the so-called anti-center direction. This field, labeled as SRa03, 
was monitored for 24.278 days starting on 1 March 2010. As a 
consequence of the DPU1 break down that took place on March 
2009, the number of targets actually photometrically monitored 
by the instrument was reduced to a maximum of 6000 stars 
only. The released telemetry is further used to oversample a 
much larger number of targets than it was initially possible, up 
to magnitude ^ 15. While not among the brightest stars of the 
field (TableQ]), CoRoT-20 benefited this new opportunity and its 
observation was performed with the regular 32 sec time sam- 
pling. It was also bright enough to allow for three-color photom- 
etry. Two trans its were detected in its light curve by the Alarm 
mode pipeline ( Sura ce et al ., 2008). The target was flagged as a 
good planetary candidate and put among the highest priorities 
for follow-up observations. 



Table 1. CoRoT-20b IDs, coordinates and magnitudes. 



C0R0T window ID 


SRa03.E2.0999 




C0R0T ID 


315239728 




USNO-B1 ID 


0902-0091920 




2MASS ID 


06305289+0013369 




GSC2.3 ID 






Coordinates 






RA (J2000) 


97.720434 




Dec (12000) 


0.22692 




Magnitudes 






Filter 


Mag 


Error 


B fl 


15.31 




V" 


14.66 




J" 


12.991 


0.023 


H" 


12.652 


0.026 


K* 


12.512 


0.027 


" from USNO-B1 


Provided by 


Exo-Dat 



The light curve of CoRoT-20b is displayed in FigQ] It shows 
a star rather quiet photometrically speaking, with no special 
prominent feature. Three transits are clearly visible with a pe- 
riod of 9.24 days and a depth slightly shallower than 1%. For the 
detailed analysis, we used the light curve reduced by the C0R0T 
calibration pipeline. It corrects for the main instrumental effects 
such as the CCD zero o ffset and gain, the back ground light and 
the spacecraft jitter (see lAuve rgne et aLl l2009l) . Portions of the 
light curve that were flagged by the pipeline as affected by parti- 
cle impacts during the South Atlantic Anomaly crossing, were 
removed and ignored in the analysis. In total, the light curve 
consists of 56869 photometric measurements. It gives a corre- 
sponding duty cycle of 88%. 



3. Follow-up observations 

A photometric time-series of the star was obtained at the Wise 
observatory on November 14, 2010 in order to check whether 
an unkn own nearby eclip sing binary could be the source of the 
transits dDeeg et alll2009l) . The detection of a transit ingress ex- 
cluded this configuration at the spatial resolution of Wise. The 
observed time of the ingress, with first contact at 24555 15.5 10± 
0.007 HJD was then used to refine the period of CoRoT-20b, to- 
wards the value quoted in Table [3] Ground-based images from 
both Wise and the DSS show that the star is rather isolated 
(Fig IS- This supports the very low contamination rate that was 
derived for the star within the C0R0T photometric mask (see 
Sec©. 
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Fig. 2. Image of the DSS showing CoRoT-20 and its environ- 
ment. The photometric mask used for C0R0T observations is 
overplotted on the target. 



Radial velocity (RV) observations started during the same 
season on December 9, 2 010. We used the HARPS spectro- 
graph dMavor et all 120031) mounted on the 3.6-m ESQ tele- 
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harps+sophie+fies Table 2. Log of radial velocity observations 




0.4 0.6 
phase 

Fig. 3. The phase-folded radial velocity measurements of 
CoRoT-20. The various symbols correspond to the different 
spectrographs used for the follow-up campaign: HARPS (black 
circle), SOPHIE (open circle) and FIES (open triangle). The best 
fit solution is over-plotted in full line 
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Fig. 4. Bisector span versus radial velocity of CoRoT-20 showing 
no correlation. 



Date 


HJD 


Vrad 


o"v rad 


Spectrograph 






kms~' 


kms~' 




2010-12-09 


55540.70645 


60.728 


0.0169 


HARPS 


2010-12-14 


55545.75533 


60.353 


0.0138 


HARPS 


2011-01-12 


55574.60927 


60.453 


0.0119 


HARPS 


2011-01-16 


55578.64264 


60.985 


0.0285 


HARPS 


2011-01-21 


55583.64035 


60.443 


0.0334 


HARPS 


2011-01-28 


55590.66159 


60.320 


0.0234 


HARPS 


2011-01-16 


55578.41934 


60.813 


0.0385 


SOPHIE 


2011-01-17 


55579.41057 


61.048 


0.0532 


SOPHIE 


2011-02-04 


55597.44200 


60.948 


0.0213 


SOPHIE 


2011-02-05 


55598.38017 


61.077 


0.0175 


SOPHIE 


2011-02-06 


55599.39134 


60.267 


0.0183 


SOPHIE 


2011-01-06 


55568.55091 


60.611 


0.032 


FIES 


2011-01-08 


55570.59973 


61.090 


0.060 


FIES 


2011-01-18 


55580.60213 


60.260 


0.036 


FIES 


2011-01-19 


55581.56506 


60.146 


0.034 


FIES 



trum of the RV standard star HD 50692 (lUdrv et all Il999h . ob- 
served with the same instrument set-up as CoRoT-20. 

The 15 radial velocities of CoRoT-20b are listed in Table [2] 
and displayed in Fig [3] They present a clear variation, in phase 
with the CoRoT ephemeris and consistent with a companion in 
the planet-mass regime with an eccentric orbit. We nevertheless 
investigated the possibility of an unresolved eclipsing binary be- 
ing the source of observed transits. To that purpose, we per- 
formed the line-bisector analysis of the CCFs (see Fig. |4]i and 
also checked that there is no dependency of the RVs variations 
with the c ross-correlation ma sks constructed for different spec- 
tral types dBouchv et al.ll2009l) . 

The Keplerian fit of the RVs was performed simultaneously 
with the transit modeling (see Section 4). All the parameters of 
the fit are listed in Table [3] 



scope (Chile) as part of t he ESQ large program 184.C-0639, 
the SOPHIE spectrograph dPerruchot et all l2008h on the 1.93- 
m telescope at the Observatoire de Haute Provence (France) 
and the FIES spectro g raph o n the Nordic Optical Telescope 
dFrandsen & Lindberg, 1999) based on the 2.56-m Nordic 
Optical telescope in La Palma (Spain) under observing program 
P42-216. 

We used the same instrument set-up as for previous CoRoT 
candidates follow-up : high resolution mode for HARPS and 
high efficiency mode for SOPHIE without acquisition of the 
simultaneous thorium-argon calibration, the sec ond fiber being 
used to monitor the Moon background light (San terne et all 
1201 ll) . For HARPS and SOPHIE, the exposure time was set to 
1 hour. We reduced data and computed RVs with the pipeline 
based on the weight ed cross-correlation function (CCF) us ing a 
numerical G2 mask (Ba ranne et al ., 1996: lPene et a Tl l2002l) . 

FIES observations were performed in high-resolution mode 
with the 1.3 arcsec fiber yielding a resolving power R » 67 000 
and a spectral coverage from 3600 to 7400 A. Three consec- 
utive exposures of 1200 sec were obtained for each observa- 
tion. Long-exposed ThAr spectra were acquired r ight before and 
after e ach science spectra set, as described in Buchhave et alJ 
(2010). Standard IRAF routines were used to reduce, combine, 
and wavelength calibrate the nightly spectra. RV measurements 
were derived cross-correlating the science spectra with the spec- 



4. System parameters 

We calculated the flux contamination from nearby stars whose 
light might fall inside the CoRoT photome tric m a sk usi ng the 
same methodology as described in iBorde et al. | (EoiOh . The 
method takes into account the photometric mask used to perform 
the on-board photometry and all the stars in the target neighbor- 
hood, including faint background stars. We found the contami- 
nation being less than 0.6% and we further neglected it. 

Three sections of the light curve, each centered on a tran- 
sit, were locally normalized by fitting a third-degree polynomial. 
Each section was a 5-hours interval before the transit ingress 
and after its egress. The detailed physical modeling of the sys- 
tem was performed by carrying out the transit modeling and 
the Keplerian fit of the radial velocity measurement s simulta- 
neously. For the transit fit we used the formalism of iGimenezI 
(20061 120091) . The fit implies twelve free parameters : the orbital 
period P, the transit epoch T tr , the transit duration Tu, the ra- 
tio of the planet to stellar radii R p /R+, the inclination i between 
the orbital plane and the plane of the sky, the Lagrangian orbital 
elements h — e sin a> and k — e cos oj, where e is the eccen- 
tricity and a) the argument of the periastron, the radial-velocity 
semi-amplitude K, the systemic velocity y re i and the two off- 
sets between SOPHIE and HARPS radial velocities and SOPHIE 
and FIES. For the transit modelin g, we used a limb-darkening 
quadratic law (IClaretl I2003L 120041) . The limb-darkening coeffi- 
cients u a and Ub were taken using the tabulated values for the 



M. Deleuil: CoRoT-20b: A very high density, high eccentricity transiting planet 



CoRoT bandpass from lSingl (120101) for the atmospheric parame- 
ters T e ff, log g and metallicity derived for the central star (see 
sect. 0>: u a = 0.4262 + 0.0168 and u b = 0.2434 + 0.0108). 
The two corresponding non-linear limb-darkening coefficients 
are u+ — u a + Uj, — 0.6696 + 0.0201 and m_ = u a - Uj, — 
0.1828 ± 0.0201. We decided to keep these limb-darkening pa- 
rameter values fixed in the transit fitting. 

The fit was p erformed using the algorithm AMOEBA 
(iPress et all 119921) . The initial values of the fitted parameters 
were changed with a Monte-Carlo method to find the global 
minimum of the x 1 - The associated 1-sigma errors were then 
estimated using a b ootstrap procedure described in details in 
iBouchv et al.l(l201 ll) . In such a procedure the limb-darkening pa- 
rameters were allowed to vary within their error bars related to 
the atmospheric parameter uncertainties. The final values of the 
fitted parameters and the subsequently derived system parame- 
ters are given in Table[3] Fig. [5]displays the best fit compared to 
the observed folded transit. 

1 .02 F : 



1 .01 : 




-0.015 -0.010 -0.005 0.000 0.005 0.010 0.015 
Orbital phase 

Fig. 5. The phase-folded transit in the phase space. The phase 
bins are 3.3 min and the error bar of each individual bin was 
calculated as the dispersion of the points inside the bin, divided 
by the square root of the number of points per bin. The best 
model is over plotted in full line. 



5. Stellar parameters 

The spectroscopic analysis was done the usual way it is carried 
out for the CoRoT planets: a master spectrum was created from 
the co-addition of spectra collected for the radial velocity mea- 
surements of the companion. We chose the HARPS spectra that 
offer the highest spectral resolution. We selected those that were 
not affected by the Moon reflected light at the time of the ob- 
servations. Each order of the selected spectra was corrected by 
the blaze, set in the barycentric rest frame and rebinned to the 
same wavelength grid with a constant step of 0.01 A. The spectra 
were then co-added order per order. Each order of the co-added 
spectrum was then carefully normalized and the overlapping or- 
ders were merged resulting in a single ID spectrum. This master 
spectrum has a S/N of 176 per element of resolution at 5760A in 
the continuum. 

A prior estimate of the atmospheric parameters T e fr, log g, 
chemical composition and v sin i was performed by fitting the 
spectrum with a library of synthetic spectra calculated using 
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Fig. 6. Abundances of the chemical elements measured with 
VWA in the HARPS co-added spectrum of CoRoT-20. The abun- 
dances refer to the solar value. White circles correspond to neu- 
tral lines, red boxes to singly ionized lines and the yellow area 
represents the mean metallicity within one sigma error bar. 



MARCS stellar atmosphere, including the Ha Balmer line. The 
rotational broadening was estimated on a selection of isolated 
spectral lines fitted by synthetic spectra convolved with vari- 
ous rotational velocities. We found v sin i = 4.5 + 1 .0 km s _1 and 
■v 'macro = 3.5 + 1 .0 km s _1 . The detailed analysis was then carried 
out using the Versatile W avelength Analysis package (VWA) 
(iBruntt et all 120041 1201 Obi) . A first set of weak and isolated lines 
of Fe i and Fe n was fitted until the derived abundances of Fe 
minimized the correlation with the equivalent width and the ex- 
citation potentials. We found : T eff = 5880 ± 90 K, log g = 4.05 
+ 0.17 and v mlc = 1.10 ± 0.1 kms 1 which corresponds to a G2- 
type dwarf. Then the abundances of other elements for which we 
could find isolated spectral lines were derived (see Fig [5]). We 
performed an independent estimate of the surface gravity from 
the pressure-sensitive lines: the Mgilb lines, the Nai D doublet 
and the Cai at 6 122 A and 6262 A. We fitted the spectrum with 
the aforementioned grid of synthetic spectra in regions centered 
on each of the spectral lines of interest. The inferred value of the 
surface gravity is log g = 4.2±0.15, a value in good agreement 
with the log g derived with VWA obtained from the agreement 
between the Fei and Fen abundances. We thus adopted log g- 
4.2 for the surface gravity. 

The mean metallicity of the star is computed as the mean of 
metals with more than 10 lines in the spectrum, such as Si, Ca, 
Ti, Fe, Ni (Fig.©. This yields a straight mean of [M/H] = 0.14+ 
0.05. The error on [M/H] due to the uncertainty on T e g , log g and 
microturbulence is 0.1 1 d ex, which we must ad d quadratically to 
get [M/H] = 0.14 ± 0.12 (IBruntt et all l20T0ah. 

We also checked for any indicators of age. We found no hint 
of stellar activity in the Ca n H and K lines. However, the Li i line 
is clearly detected at 6708 A (see Fig. [7]). We measured an equiv- 
alent width of W ea = 44 mA and determined a lithium abundance 
of 2.97. Following lSestito & Randichl (12005). this leads to an age 
in the range 100 Myr up to 1 Gyr, depending on the star's initial 
rotation velocity. 

The modeling of the star in the HR diagram was carried out 

1 /3 

in the (T e ff,MJ //?*) plane taking the host-star's metallicity into 
account. It resulted in the final estimates of the star's fundamen- 
tal parameters given in Table[3j M* = 1 . 14 ± 0.08 M e , R+ = 1.02 
± 0.05 R Q . The inferred surface gravity is log g = 4.47 + 0.11, 
in agreement within the errors with the spectroscopic result. The 
evolutionary status points to a young star likely in the last stages 
of the pre-MS phase. We found the most likely isochrone age 
being 100^2° Myr, a result in good agreement with the Li abun- 
dance. 
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Wavelength [A] 

Fig. 7. CoRoT-20 spectrum in a spectral region around the Lii 
lines at 6708A. 



We calculated the distance of the star. We used the parame- 
ters of the star we derived and its 2-MASS magnitudes to estimate 
the reddening. We found a color excess E( J - K) — 0. 1 8 mag and 
the absorption Ay - 1 ,04±0.5 mag using the extinction law from 
ISchlegel et ail (119981) . This yields a distance of 1.23 + 0.12 kpc, 
consistent with the strong interstellar absorption observed in the 
Nai(Dl)and (D2) lines. 



6. CoRoT-20 system properties 

Compared to the sample of known transiting planets, CoRoT- 
20b is unusual in many respect. With an orbital period of 9.24 
days it joins the group of transiting planets with periods out- 
side the pile-up at 3 days. It is the fifth planet discovered by 
CoRoT in this period domain which currently accounts for 25 
planets (see http://exoplanet.eu/), 9 out of the s e bel onging to 
multi -planet systems : K epler-9 (Holman et al., 2010J), Kepler - 
10 dFressin et all 120111) and Kepler- 11 (llissauer et all 120111) . 
However all these Kepler-planets have a mass which is less than 
~ 0.3 Mj up and could not be directly compared to the giant planet 
population. Excluding these planets in multiple systems, for the 
17 remaining objects of the sample that do not have a detected 
companion, 8, that is 47%, have a significant eccentric orbit with 
e in the range 0.15 to 0.9. 

Planets with highly eccentric orbit appear to be found pref- 
erentially among the high-mass and/or long period planet popu- 
lation. With a mass of 4. 13 Mj up which places it at the border of 
the gap in mass between the regular hot-Jupiter population and 
the very massive planet one, CoRoT-20b is consistent this trend. 
In the mass-period diagram they ar e clearly separated from the 
lighter planets with circular orbits dPont et all 1201 ll) . This di- 
chotomy and in particular the lack of massive close-in planets at 
circular orbit suggest that tidal evolution should play an impor- 
tant role in the fate of the planet population. 

6.1. Tidal evolution 



Following iLevrard et alJ d2009l) approach we checked the sta- 
bility of CoRoT-20b to tidal dissipation. The authors calculated 
the ratio between the total angular momentum of a given system 
L tot and the critical angul ar momentu m L cr „ for some transit- 
ing systems. According to Hut (1980), tidal equilibrium states 
exist when the total angular momentum is larger than this crit- 
ical value L cr j,. However, this equilibrium state could be stable 
or unstable, depending whether the orbital angular momentum 
L„,-b is more than three time the total spin angular momentum 
L sp in, or not. lLevrard et al.l (120091) demonstrated that for none of 
the systems but HAT-P-2b the stable tidal equilibrium state, that 
corresponds to L tot /L cr i t > 1, exists. Further the fate of these 




Fig. 8. Tidal evolution of the orbital semi-major axis and eccen- 
tricity. The figure is displayed on a time interval larger than the 
expected lifetime of the star to show the triple synchronization 
characteristic of a Darwin-stable system. 

close-in planets is ultimately a collision with their host-star. The 
study has been recently reexamined and extend ed to more than 
60 transiting systems by Matsumur a et alJ 12010) who achieved a 
similar conclusion, showing that the vast majority of these close- 
in planets will spiral-in to their host star an d will be destroyed by 
tides. Using equations (1) and (2) given bv lLevrard et al.l (2009) 
that neglect any effect of a possible magnetized stellar wind, we 
found for CoRoT-20b: 



Ltoi I L c 



1.057 and L 



spin* I ^orb 



0.0458 



It shows that, within the current observational uncertainties, the 
planet has a tidal equilibrium state. It is worth noticing that our 
approach also assumes the stellar obliquity is small. The later is 
poorly constrained as the star's rotation period could not be de- 
rived from the light curve. We simply assumed that the rotation 
axis is perpendicular to the line of sight and derived the star's 
rotation period from the v sin i (Table [3), a regular method for 
transiting systems. This gives a rotational period of the star of 
1 1.5 ± 3.1 days, that is of the same order than the planet's orbital 
period. In the case of CoRoT-20b, L sp i ni ,/L or t, < 1/3 and most of 
the angular momentum o f the system is in the orbit. According to 
iMatsumura et a l. (2010), CoRoT-20b belongs to the very small 
subgroup of Darwin-stable systems that evolve toward a stable 
tidal equilibrium state where migration will stop. 

From the Roche-limit separation, the planet th us lies well 
beyon d two times the Roche limit distance. Using dFaber et all 
120051) : 

a R = (R p /0A62)(MJM p ) l/3 

we found that the Roche limit a« of the system is 0.0057 
AU. This further supports the migration scen ario over the 
scattering/Kozai-cycle scenario as proposed by iFord & R asio 
d2006h . 

We performed a complete calculation of the tidal evolution 
of the system for med by the star and t he pla net assuming a 
linear tidal model dMignardl 119791; iHuil 119811) . The main dif- 
ficulty here is to choose the values of the dissipation in the star 
and in the planet. For the main semi-diurnal tides of the star, 
we have adopted the value Q' = 10 7 as foun d for hot Jupiters 
dHansenl 12010b iBenftez- Llambav et all 1201 ll) . Because of the 
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close values of the orbital period and the rotation of the star, 
the components of the tides raised on the star by the planet re- 
lated to the orbital eccentricity are also equally important, but 
the values of the current dissipation obtained with them are of 
the same order of magnitude. For the planet, we have derived 
one value on the basis of the actually deter mined Q' of Jupiter 
(Q = 1.36 x 10 5 , see lLainev eFaD (120091)1 . We first note that 
standard linear tidal theories (see lHutl fl981. eqn 45) allow us to 
determine the current rotation period (stationary) of the planet 
independently of the dissipation. We obtain 2.64 + 0.13 days. 
To transform Qt into the planet's Q' p , we have to take into 
account: (i) Q'„ sc ales with the semi-diurnal tide period (see 
IFerraz-Mello et alll2008l:lMatsumura et aUl2010l) ; (i i) Q' p s cales 
with R~ 5 (see lEggleton et all 119981; lOgilvie & LinT 120041) . We 
thus obtain Q' p = 2.2 x 10 6 . Fig. [8] shows the variation of the 
semi-major axis and the evolution of the eccentricity. One can 
notice that with the adopted dissipation values, while the eccen- 
tricity tends toward zero, the circularization will not be achieved 
within the lifetime of the system. 

Fig. [9] shows the evolution of the periods. The planet ro- 
tation is currently in a stationary super-synchronous state, that 
is the planet rotation is faster than its orbital motion. Its period 
increases as the eccentricity decreases and almost synchroniza- 
tion is reached when the eccentricity becomes very small. The 
star rotation period is currently decreasing; it will equal the syn- 
chronous value at some time 4 Gyr from now. However, it will 
continue to decrease up to reach the triple synchronous station- 
ary state. The triple synchronization, however, does not seem to 
be reached within the lifetime of the star. 

It is worth underlining that the actual Q' values are not 
known and the values we used are only estimates founded on 
previous studies. Therefore the exact timescale of the tidal pro- 
cesses is uncertain. Furthermore, by extracting angular momen- 
tum from the system, stellar magnetic braking may prevent the 
planet from reaching a triple synchronous state and ultimately 
jeopardize its survival dBouchv et al.Ll20lTl) . Indeed, simulations 
in which magnetic braking was active duri ng the whole system 
lifetime, following the model proposed by (iBouvier et all 09971) 
and using the same tide parameters as in the examples given 
above, show that the planet is falling below the Roche limit in 
about 6 Gyr. This result is critically dependent on the adopted 



parameters and further would required a detailed study that is 
well beyond the scope of the present paper. 

We also investigate the consequences of the circularization 
of the planet orbit which is in the phase of fast circularization, on 
the transits occurrence. Assuming there is no other close massive 
perturber in the system, then two effects are causing TTVs: the 
decrease of the orbital semi-major axis and the circularization of 
the orbit. Concerning the orbital semi-major axis, the time-scale 
of its variation, a, is -0.95 10 -5 1/Myr presently (see Fig [9j>. As 
a consequence, a continuous period variation of P/P « -410" 12 
per cycle is expected. As well-known, this linear period variation 
will cause a parabolic O—C curve, and in 100 years from now the 
O - C value will be only -25 sec onds. This i s slightly over the 
3cr ob servation limit by CoRoT (Beanl 120091; ICsizmadia et all 
2010). Assuming that the transit timing precision can be forced 
down to 5 sec in the future, this O - C value will be reached in 
45 years from now. 

Turning to the evolution of the eccentricity during the cir- 
cularization process, it has two consequences. First the occur- 
rence of the secondary eclipse will change. The displacement 
D of the s econ dary from phase 0.5 is given by (eqn 1 and 2 
iBorkovitsl 12004 e.g.). The previous results of the tidal evo- 
lution calculations indicated e = -4.5 10~ 5 1/Myr and P = 
-1.5 10~ 3 days/Myr. Assuming a constant to, we have that 
D = -37.56 10~ 5 days/Myr or D = -9.53 10~ 12 days/cycle. 
This variation is of the same order as the previous one caused by 
the decreasing semi-major axis, so it would be observable within 
a century, too. 

For the second effect, that is the circularization of the or- 
bit, one can also consider the occurrence of a small precession 
of the orbit. This effect is hardly observable, but interesting on 
the theoretical side, since the transit occurs at the true anomaly 
v = 90° - to where a> is the argument of periastrion. The later is 
also subject to variations because of theory of general relativity 
but also because the tidal effects force the apsidal line to rotate. 
However, this variation has a different time-scale. We thus do 
not take this into account here, even if tidal forces also cause a 
small precession of the orbit showing that co is not zero. So if e 
decreases due to circularization, and even if to is constant, then 
at the epoch of transit the eccentric anomaly will increase and 
hence the mean anomaly at transit will occur later. However, a 
first estimation shows that this effect may be negligible in a ten 
year timescale. 

6.2. Internal structure 

CoRoT-20b is a massive hot-Jupiter with a mass of 4.24 Mj up 
a radius of 0.84 Mj up , and an inferred density p = 8.87 ± 
1.1 gem" 3 . A few giant planets are already reported with sim- 
ilar density or even higher : CoRoT-14b p = 7.3 + 1.5gcm~ 3 
dTinglev et all 1201 ll). WA SP- 18b p = 8.8 ± 0.9 dHellier et all 
120091: ISouthworthi 120 1 Oh or HAT-P-20b p = 13.78 + 1.5 
g cm 3 dBakos et all 1201 Oh for example. While the mass of these 
planets spans a large range, from ~ 4 up to more than 9Mj llp , 
their radius is close to lRj up . Given CoRoT-20b\ large plane- 
tary mass, its small size is surprising. Among these high den- 
sity giants planets, only HAT-P-20b has a comparable size, i.e. 
0.867 ± 0.033 R Sup . CoRoT-20b, as HAT-P-20b, is thus expected 
to contain large amounts of heavy elements in its interior. 

To investigate the internal structure of CoRoT-20b, 
we computed p la netary evolution models with CEPAM 



-H ) 1 r 

(Guillot & Morel, 1995), following the description in 



Guillot & HaveT(l201 ll) . and lHavel et all (1201 ll) for a planet of 
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a total mass 4.24 Mj up . We derived a time-averaged equilibrium 
temperature of the planet to be T eq = 1002 ± 24 K. The 
results for T eq = 1000K are shown in Fig. [10] in terms of the 
planetary size as a function of the system age. The coloured 
regions (green, blue, yellow) indicate the constraints derived 
from the CESAM stellar evolution models ( M orel & Lebretonl 
2008) at 1, 2, and 3cr level, respectively. For preferred ages 
between 100 Ma and 1 Ga, we find that CoRoT-20b should 
contain between 680 and 1040 M e of heavy-elements in its 
interior (i.e. between 50 and 77% of the total planetary mass), 
at lc r level, about twice the amount needed for HAT-P-20b 
(see iLeconte et al.L 1201 lb . While this result is qualitatively 
in line with the observed correlation betw een star metallicity 
and heavy elements in the planet (e.g. iGuillot et all [2006; 
Mille r & Fortnevt 1201 1[ and references therein), the derived 
amounts are extremely surprising. They would imply that all the 
heavy elements of a putative gaseous protoplanetary disk of 0. 1 
to O.15M were filtered out to form CoRoT-20b, and then that 
an extremely small fraction of hydrogen and helium in that disk 
was accreted by the pl anet. This i s at odds with todays accretio n 
formation models (e.g. llda & LinLl2004tlMordasini et al., 2009). 
Disk instability with di fferentiation and partial tidal stripping 
dBolev & Durisenll2010) is proposed as an alternative formation 
pathway. According to Bolev et ail (1201 ll) . this could account 
for giant planets with massive core such as HAT-P-20b. In its 
current orbit, CoRoT-20b doesn't enter the Roche limit and no 
tidal stripping is acting on yet but this scenario would deserve 
further investigations. 

We investigated the possibility that changes in the atmo- 
spheric model would yield more "reasonable" values for the 
planetary enrichment . As can be seen from a similar study in the 
brown dwarf regime (Bu rrows et al.Ll20Tlb . the consequences of 
modified atmospheric properties are limited for objects with the 
mass of CoRoT-20b (i.e. standard radii for objects of this mass 
range from 1.05 to 1.20 Rj up ). By artificially lowering the in- 
frared atmospheric opacity by a factor 1000 (not shown), we 
were able to decrease the lcr upper limit to the core mass from 
650 to 390 M ffi , a small change compared to huge and unphysical 
decrease in the opacity. 

On the other hand, one strong assumption in our study is that 
heavy elements are embedded into a central core. When rela- 
tively small amounts of heavy elements are considered, it is not 
very important whether they are considered as bein g part of a 
core or mixed in the envelope (e.g.llkoma et al.. 2006). However, 
as shown by [Baraffe et al. (2008), when 0.5 Mj up of ices are 
mixed in the envelope of a 1 Mj up planet, its radius is smaller 
by ~ 0.1 Rj np than when one considers that these elements are 
part of a central core. It is thus very likely that the mass of heavy 
elements required to explain the radius of CoRoT-20b is high 
but si gnificantly smaller th an considered here. Estimates based 
on the lBaraffe et al .1 (120081) calculations indicate that if mixed in 
the envelope, a mass of heavy elements 2 to 3 times smaller than 
estimated in Fig. [10]would explain the observed planetary size. 
This would alleviate the problem of the formation of the planet, 
although it would still require relatively extreme/unlikely sce- 
narios. 

7. Summary 

In this article we presented the discovery of CoRoT-20b. The 
object belongs to the population of massive planets with or- 
bital semi major axes below 0.1 AU, a domain of orbital peri- 
ods where low and high eccentricity systems co-exist in a nar- 
row range of orbital period. We examined the tidal stability of 
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Fig. 10. Evolution of the size of CoRoT-20b (in Jupiter units) as 
a function of age (in billion years), compared to constraints in- 
ferred from CoRoT photometry, spectroscopy, radial velocime- 
try and stellar evolution models. Green, blue and yellow regions 
correspond to the planetary radii and ages that result from stellar 
evolution models matching the inferred p* - T e ff - [Fe/H] un- 
certainty ellipse within lcr, 2cr and 3<x, respectively. Planetary 
evolution models for a planet with a solar-composition envelope 
over a central dense core of variable mass (0, 400, 800, and 1000 
M e as labelled) are shown as dashed lines. These models also 
assume that 1 % of the incoming stellar irradiation is dissipated 
deep into the interior of the planet. 



CoRoT-20 and found that, within the observational uncertain- 
ties, it belongs the relatively small population of transiting plan- 
ets that are considered as "Darwin-stable", i.e. systems for which 
in the absence of processes extracting angular momentum from 
the system (i.e. stellar winds), the planet would never fall onto 
the central star which would instead be spun-up towards triple- 
synchronization (equality of the orbital, planetary and stellar 
spin periods). Other cases in the planeta ry domain are CoRoTSb , 
CoRoT-6b, WASP-7 and HD 80606 (Mat sumuraet all [2010). 
Measuring the stellar obliquity through the Rossiter-McLaughlin 
effect would provide additional constraints to better constrain 
its tidal evolution and further understand its orbital evolution. 
The expected semi-amplitude of the radial velocity anomaly of 
the Rossiter-McLaughlin effect is estimated to be 22 + 5 ms~', 
unfortunately quite difficult to detect with the present spectro- 
graphs used for the keplerian orbit determination. Another point 
would be to assess the presence of any additional companions in 
the sy stem by long-term radial velocity monitoring as, accord- 
ing to Matsu mura et al.l d2010h . the formation path of close-in 
planet should be different for single-planet system and multi- 
planet ones. CoRoT-20 system appears thus an interesting bench 
test case to study the tidal orbital and rotational evolution of the 
close-in population. 

The second interesting peculiarity that distinguishes CoRoT- 
20b from the regular giant planet population is its small observed 
radius. According to planetary evolution models, the interior of 
this compact planet should contain a very high amount of heavy 
elements, with a central dense core whose mass would be in the 
range between 680 and 1040 M & . Although mixing heavy ele- 
ments in the envelopes rather than confining them to a central 
core can lead to substantially smaller values (by a factor esti- 
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mated to be ~ 2 - 3), the origin of such a huge amount of heavy 
elements is difficult to explain within the framework provided 
by the current pla netary formation models. With HAT-P-20b 
(Bak os et al.Ll2Qldh . it is the second example of extremely metal 
rich interior that challenges planetary interior models. However, 
the two planetary systems differ in many aspects: HAT-P-20 
mass is nearly twice that of CoRoT-20b; it orbits on a nearly cir- 
cular orbit a K3 metal rich star, while CoRoT-20 is a solar-type, 
slightly metal enriched star. In addition, HAT-P-20 seems to be 
physically associated to another stellar companion, while up to 
now, CoRoT-20 has none, detected or suspected. It thus appears 
difficult from this two exceptions to derive any trend that would 
provide hints on the origin these challenging and intriguing bod- 
ies. 
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Table 3. Planet and star parameters. 



Ephemeris 



Planet orbital period P [days] 
Primary transit epoch T tr [BJD-2400000] 
Primary transit duration T\\ [d] 
Secondary transit epoch T s [B1D-2400000] 

System parameters 



9.24285 ± 0.00030 2 
55 266.0001 ± 0.0014 3 
0.0927 ± 0.0019 
55 272.46 ±0.13 4 



Periastron epoch T peri [BJD-2400000] 


55265.79074 


5 


esinoi 


0.468 ± 0.017 




ecosw 


0.312 ± 0.022 


6 


Orbital eccentricity e 


0.562 ± 0.013 




Argument of periastron oj [deg] 


56.3 + ^ 


7 


Radial velocity semi-amplitude K [ms~'] 


454 ± 9 




Systemic velocity y re ; [ km s~ 1 ] 


60.623 ± 0.006 


8 


HARPS-SOPHIE off set velocity V r j [ms" 1 ] 


93 ± 11 




SOPHIE-FIES offset velocity V r2 [ms^ 1 ] 


163 ± 20 


9 


O-C residuals [ m s~'] 


27 


10 


Radius ratio Rp/R* 


0.0842 ± 0.0017 




Impact parameter" b 


0.26 ± 0.08 


1 1 

12 


Scaled semi-major axis a/R+ h 


18 95+ - 63 


M^/R* [solar units] 


1 

-0.039 


13 


Stellar density p* [g cm~ 3 ] 


1 51 +015 




Inclination i [deg] 


88.21+ 0.53 


14 


Spectroscopic parameters 







Effective temperature T e ff[K] 


5880 ± 90 


15 


Surface gravity log g [dex] 


4.20 ±0.15 




Metallicity [Fe/H] [dex] 


0.14 ±0.12 


16 


Stellar rotational velocity v sin i [ km s~ 1 ] 


4.5 ± 1.0 


Spectral type 


G2 V 


17 


Stellar and planetary physical parameters from 


combined analysis 


Star mass [M Q ] 


1.14 ±0.08 


18 


Star radius [R Q ] 


1.02 ± 0.05 




Distance of the system [kpc] 


1.23 ±120 


19 


Age of the star t [Myr] 


100™ 


Orbital semi-major axis a [AU] 


0.0902 ±0.0021 


20 


Orbital distance at periastron a p „ [AU] 


0.0392 ± 0.0017 


Orbital distance at apastron a apo [AU] 


0.1409 ±0.0037 


2] 


Planet mass M p [MjY 


4.24 ± 0.23 


Planet radius R p [Rj] c 


0.84 ± 0.04 




Planet density p p [g cirT 3 ] 


8.87 ± 1.10 




Equilibrium temperature d T eq [K] 


1002 ± 24 




Equilibrium temperature at periastron ' Tj,'" [K] 


1444!" 




Equilibrium temperature at apastron'' [K] 


764 ±18 
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a / K * ~ l_ e 2 

lGimenezLl2009l) . 

b 



^l-cos-(v') +oj- ^)-sin- i 



, where v\ is the true anomaly measured from the periastron passage at the first contact (see 



qgos / 
R, 



1+e-sin uj 

Radius and mass of Jupiter taken as 71492 km and 1.8986xl0 30 g, respectively. 
d zero albedo equilibrium temperature for an isotropic planetary emission. 
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